Asymptotic Euler-Maclaurin formula over lattice polytopes 



Tatsuya Tate* 
Graduate School of Mathematics 
Nagoya University 
Furo-cho, Chikusa-ku, 
Nagoya, 464-8602, Japan 
Emaik tate@math.nagoya-u.ac.jp 



Abstract 

Formulas for the Riemann sums over lattice polytopes determined by the lattice points 
in the polytopes are often called Euler-Maclaurin formulas. An asymptotic Euler-Maclaurin 
formula, by which we mean an asymptotic expansion formula for Riemann sums over lattice 
polytopes, was first obtained by Guillemin-Sternberg [GSj . Then, the problem is to find a 
concrete formula for each term of the expansion. In this paper, an asymptotic Euler-Maclaurin 
formula of the Riemann sums over general lattice polytopes is given. The formula given here is 
an asymptotic form of the so-called local Euler-Maclaurin formula of Berline-Vergne |BeVj . For 
Delzant polytopes, our proof given here is independent of the local Euler-Maclaurin formula. 
Furthermore, a concrete description of differential operators which appear in each term of the 
asymptotic expansion for Delzant lattice polytopes is given. By using this description, when 
the polytopes are Delzant lattice, a concrete formula for each term of the expansion in two 
dimension and a formula for the third term of the expansion in arbitrary dimension are given. 

Introduction 

In this paper, we consider asymptotic behavior of the Riemann sums over lattice polytopes, 

R^(p-^^y-=j;^ E ^W^)' (0-1) 

76(AfP)nZ'" 

where P is a lattice polytope in M"^, which means that each vertex has integer coordinates, and ip is 
a smooth function on P. Formulas for R]\f(P; ip), which are often called Euler-Maclaurin formulas, 
are extensively investigated in combinatorics and geometry of toric varieties. If we take = 1, the 
Riemann sum Rn{P', 1) is reduced to the so-called Ehrhart polynomial 

Ep{N) := tt(iVP) n = N'^''"^^^Rn{P; 1), 

which is closely related to the Todd class of a toric variety corresponding to the polytope P. 
In this context, geometry of toric varieties is a suitable and powerful tool to analyze the function 
Ep{N). Indeed, as in |Fuj . one can show that Ep{N) is a polynomial in by using the Hirzebruch- 
Riemann-Roch theorem. The problems concerning (exact) Euler-Maclaurin formulas and Ehrhart 
polynomials are investigated by various authors, for example |CS| . |BeV| . [BrV] . |KP] . See [KSW2] 
and references therein for various results on these topics. 

Before explaining some of the results closely related to the present paper, we state one of our 
theorems. 
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Theorem 1 Let P he a lattice polytope in M™. For each face f of P and non-negative integer 
n with dim(/) > dim(P) — n, there exists a homogeneous differential operator Dn{P; f) of order 
n — dim(P) + dim(/) with rational constant coefficients which involves derivatives only in directions 
orthogonal to the face f such that for each smooth function ip on P, we have the following asymptotic 
Euler-Maclaurin formula: 

i?^(P;</.) ~ J^iV-" Yl [Dn{P;f)v (iV^oo), (0.2) 

where F{P) denotes the set of faces of P. The integration in the right hand side is performed with 
respect to the measure on the affine hull ( / ) of f which is the parallel translation of the Lebesgue 
measure on the subspace L{f) parallel to {f) defined by the lattice L{f) n A. 

In this section, we explain some of the previous works on the Euler-Maclaurin formula closely 
related to Theorem [1] and mention other results obtained in the present paper. 

An exact Euler-Maclaurin formula for Delzant polytopes was originally obtained by Khovanskii- 
Pukhlikov |KPj . and Brion-Vergne [BrV^ generalized it to simple polytopes without using the theory 
of toric varieties. One of their results can be stated as (assuming that P is a Delzant polytope) 



Ri^{P-^) =Todd{P;d/Ndh) [ <f{x)dx 

JPh 



(0.3) 

h=0 



where (p is a polynomial, h = {hi, . . . , hd) G is a small parameter with d the number of faces of 
P of codimension one, Todd(2:) = ^^^-z is an analytic function around the origin, called the Todd 
function, 

d 

Todd{P; d/Ndh) = JjTodd(9/A^9/ii) 

1=1 

is a differential operator (of infinite order), and when the polytope P is given by P = {x ; ( Uj, x ) > 
Ci,i = then P^ = {x ; {ui,x) > Ci — hi,i = l,...,d}. Note that Brion-Vergne |BrV) 

obtained the same formula for simple polytopes with a modification of the differential operator 
Todd{P; d/Ndh). 

In |BeVj ■ Berline-Vergne obtained an effective formula for R]\f{P; ip) (still p) being assumed to be 
polynomial), which they call a local Euler-Maclaurin formula. This formula is of the form (setting 
N = 1 for simplicity) 

Ri{P-p,) = Y, I D{P,f)^, (0.4) 

where the sum runs over all faces / of P, D{P, f) is a differential operator (of infinite order) with 
rational constant coefficients on which involves derivatives only in directions perpendicular to 
the face /. One of remarkable points is that the formula (10. 4p of Berline-Vergne holds for any 
rational polytopes, which means that each vertex of the polytope has rational coordinates. They 
constructed a meromorphic function fi{a) for any affine rational polyhedral cone a and use a sort of 
inclusion-exclusion property (which is called a valuation property) of // to show that it is analytic 
near the origin, and they define the symbol of the operator D{P, f) by using /i. 

The operators Dn{P; f) in our formula (j0.2p is, by definition, the homogeneous parts of the 
operator D{P, f) in (j0.4p . Thus, one can think the formula (j0.2p as an asymptotic form of the local 
Euler-Maclaurin formula (10. 4p due to Berline-Vergne. As we point out in Subsection II. 3^ one can 
deduce ()0.2p by using one of results in jBeV) directly and formally. However, the method mentioned 
in Subsection 11.31 is formal, and we use a different method to prove Theorem [TJ Moreover, any 
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transparent formula for the homogeneous parts of D(P, f) is, in general, not known. We will 
see that, when P is a Delzant lattice polytope, the operators D„,(P;/) can be, to some extent, 
expressible concretely (Definition 13.61 Theorem I3.9p . Note that our formula ()0.2p is valid for any 
smooth function ip on P. Our construction of the operator Dn{P', f) makes us to obtain concrete 
formula for Delzant lattice polytopes in two dimension (Corollary 15.41) . A part of our construction 
of these operators Dn{P]f) uses an induction procedure, and they are still complicated. This 
complication comes from the "angles" at each face of the polytopes, and hence it would be rather 
natural. The complication involving the "angles" is embodied in an integration by parts procedure. 

In this paper, by the name asymptotic Euler-Maclaurin formula, we mean formulas of asymp- 
totic expansion of the Riemann sum R]\j{P]Lp). In one dimension (m = 1 and P = [0,1]), the 
following asymptotic Euler-Maclaurin formula is well know. 



i-^V'(A:/iV)=ii^([0,l];(^) 



^T^ (0.5) 

-2n 



^ / v.(x)dx+^(y(l)-^(0)) + ^ ^ ^^^"J/" (y(^"-^)(l)-^(^"-^)(0)) AT- 

•^0 n>l ^ ^' 

where (f is any smooth function on [0, 1], and 6„ are the coefficients of the Taylor expansion of the 
Todd function: 



oo , 

Or. 



Todd(-z) = 



n=0 

\n—li 



and Bn = ( — 1)"~ &2n {n > 1) are the Bernoulli numbers. 

A higher dimensional analogue of (|0.5p was given by Guillemin-Sternberg ( |GSj ). Namely, 
Guillemin-Sternberg obtained the asymptotic Euler-Maclaurin formula of the form (assuming that 
P is Delzant) 



RNiP;^) ^Todd{P;d/dNh) [ (p{x)dx 

Jpi, 



(0.6) 

h=0 

This formula also holds true for simple lattice polytopes under a modification. Note that this 
formula is, at least its appearance, similar to the Brion-Vergne formula (j0.3p . The proof of (j0.6p 
given in |GS| is different from the proof of (10. 3p given in [BrVj , and it does not use geometry of toric 
varieties. There are some applications of the above formula for spectral analysis on toric Kahler 
manifolds. In fact, in [GW| . the asymptotic Euler-Maclaurin formula obtained in [GS], combined 
with an asymptotic expansion of 'twisted Mellin transform' studied in |Wj . is applied to analyze a 
spectral measure on C™ which is, in a GIT setting, related to the pair {X, L) where X is a toric 
manifold corresponding to a Delzant polytope and L is a Hermitian line bundle on X. (See also 
[Ch] where the same spectral measure as in |GW] is discussed.) 

One more asymptotic Euler-Maclaurin formula was brought to us by Zelditch [Zj. The formula 
obtained in [Z] is stated as 

Rn{P;v)^ [ Vdx + ^ [ (/j(x) dcT + V iV-" / £n{PMx)dx, (0.7) 

Jp m Jqp ^ Jp 

where P is a Delzant lattice polytope, £n{P) is a differential operator (of finite order), and da is the 
Leray measure on the boundary dP. In [Z], Zelditch introduced the notion of Bergman-Bernstein 
measures (this name is taken from [T]) and obtained its asymptotic expansion. Then, integration 
(over the toric Kahler manifold corresponding to the Delzant polytope P) of the asymptotic expan- 
sion yields the formula (|0.7p . In [Z], the formula (|0.7p is called a 'metric expansion' to distinguish it 
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from the Euler-Maclaurin formula of the form (jU.6|) . since the differential operators £n{P) depend 
on the choice of a Hermitian metric on a line bundle over the toric manifold. But, the Riemann 
sum itself does not depend on such a metric. A point is that such a metric dependence would be 
disappeared after an integration by parts. Indeed, in |Z], the second term is computed by using an 
integration by parts identity due to Donaldson [D] . 

As is mentioned in [Z], comparison of asymptotic Euler-Maclaurin formula and the metric 
expansion of the form (10. 7p will give some further identities in the lower order terms. One of our 
motivation is to give another asymptotic Euler-Maclaurin formula which is computable to some 
extent. Indeed, we have a concrete formula for the third term of the expansion when the polytope 
is Delzant. See Corollarv 15.61 in Subsection 15.31 Thus, if one can compute the differential operator 
£2{P) in ()0.7p in terms of curvatures, then one will obtain an integration by parts identity in the 
third term in ()0.7p . which might be useful to geometry of toric manifolds. 

An idea of proof of Theorem [1] is to reduce the problem to that for unimodular cones, which 
are cones generated by a part of an integral basis, by using a subdivision of a rational cones 
into unimodular cones (see [Fuj . Section 2.6) and a canonical decomposition of the characteristic 
functions of polytopes (see the equations (|5.3p . (|5.6p in Section [5]). The asymptotic Euler-Maclaurin 
formula of Riemann sums over unimodular cones can be deduced by a method in [GS] (see also jAMj , 
|KSW1| . |KSW2] ). However, we deduce it here by a quite different method. This method is rather 
similar to the Bergman-Bernstein approach in [Z]. But, we work on unimodular cones instead of 
polytopes themselves. Thus, we use the Szasz measures introduced in Section [2] instead of Bernstein 
or Bergman-Bernstein measures discussed in [Zj or [T]. More concretely, an asymptotic property 
of the Szasz functions is used to show Proposition 13.11 in Section [Sj which is an asymptotic Euler- 
Maclaurin formula for unimodular cones. Proposition 13.11 can be deduced directly from Theorem 
3.2 in [GS], and one can consider that the Proposition 13.11 is a starting point for the subsequent 
sections. Thus, one might be able to perform similar computations in sections after Section [3] at 
least for simple polytopes, by using Theorem 3.3 in [GSj instead of Proposition 13.11 However, the 
asymptotic behavior of Szasz functions would be a general interest in its own right. Furthermore, 
there would be a possibility of using a version of Szasz functions to get asymptotics of the Riemann 
sum over general rational cones without using a subdivision of cones into unimodular cones, if one 
could resolve a problem on 'rare events' along with the lines in [T]. (See also Remark after the proof 
of Theorem l5.11 on this point.) In one dimension, we compute explicitly each term of the expansion 
for twisted Riemann sum by using this approach. This computation uses the twisted version of the 
Szasz function, and it shows that coefficients in the Taylor expansion of the 'twisted' Todd function 
can be represented by the Stirling numbers of the second kind (in particular, the equation p.lOp ). 
which is a generalization of a well-known formula among Bernoulli numbers, Catalan numbers and 
the Stirling numbers of the second kind (see ()2.17p or [GKP] ). Thus, this approach might have 
some advantages also in higher dimension. These are the reasons why we use the approach with 
the Szasz functions in this paper. 

We here mention that an asymptotic expansion of the Szasz function was first obtained in [Fe| . 
In [Fej , Feng also obtained an asymptotic formula of the Riemann sum over the positive orthant 
in the same strategy as ours. However, concrete formulas for each term of the asymptotic expansion 
are not discussed fully in [Fe]. We give an explicit formula for each term of the expansion of the 
Szasz function in Section [2j (The main purpose in |Fej was to give a non-compact analogue of 
Bergman-Bernstein approximation in [Zj. Indeed the Szasz function, defined in Section [2] in the 
present paper, is closely related to the Bergman kernel for the Bargmann-Fock space as explained 
in [E]-) 

We close Introduction with some comments on the organization of this paper. We collect some 
of the notation used in this paper in Subsection 11.11 and then, we review and define the Berline- 
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Vergne operators Dn{P] f) in Subsection 11.21 As we mentioned above, a heuristic argument to find 
a formula (j0.2p is given in Subsection ll.31 In Subsection [T31 we prove a uniqueness theorem on the 
expression of each term of the asymptotic expansion of the form ()0.2p (Theorem II. 2p . In Section 
[21 we study asymptotic behavior of Szasz functions. Some computations for the twisted Riemann 
sum in one dimension is given in Subsection 12.11 In Subsections 12.21 12.31 we define and study Szasz 
functions and their asymptotic behavior. Section [3] is devoted to the study of asymptotic behavior 
of the Riemann sums over unimodular cones. First, we prove an asymptotic expansion formula 
(Proposition 13. ip by using the asymptotic property of the Szasz functions studied in Section [2j 
Asymptotic formula obtained in Proposition 13.11 uses differential operators in direction transversal 
to each face of the unimodular cone. Then, one can perform further integration by parts. This 
is done in Subsection 13.21 In Subsection 13.31 we define differential operators obtained by the 
integration by parts procedure discussed in Subsection 13.21 which is used to renormalize each term 
of the expansion in Proposition 13. 1[ The fact that the operators so defined coincide with the 
Berline- Vergne operators is proved also in this subsection (Theorem 13. 9p . In Section HI we prove 
the asymptotic Euler-Maclaurin formula for general pointed rational cones by using the Berline- 
Vergne operators and the subdivision of pointed rational cones into a finite number of unimodular 
cones. Finally, in Section O we prove our main Theorem [H which is reformulated in Theorem l5.lt 
and a uniqueness result (Theorem 15. 3p . and give some explicit computation. 

Acknowledgment. The author would like to thank to Dr. Micheal Stolz who informed him 
about the work of O. Szasz [S]. He would also like to thank to Prof. Steve Zelditch for his helpful 
comments on the earlier version of the paper. 

1 Berline- Vergne operators and heuristic argument 

In this section, we review the symbol of differential operators defined in |BeV] . Then, we give a 
heuristic argument to obtain an asymptotic Euler-Maclaurin formula of the form (|0.2p . Further- 
more, we deduce a uniqueness theorem on expression of coefficients in asymptotic Euler-Maclaurin 
formula of the form (10. 2p . 

1.1 Notation 

Let X be a finite dimensional vector space over M, and let A be a lattice in X. Such a pair (X, A) is 
called a rational vector space. The dual space X* of a rational space {X, A) is a rational space with 
the dual lattice A* of A. A point x £ X is said to be rational if G A for some q € Z\{0}. The set 
of rational points in X is denoted by Xq. A basis of A over Z is called an integral basis of A. For 
each rational vector space {X, A) , we fix a Lebesgue measure on X normalized so that the measure 
of the fundamental domain of the action of A on X has measure 1. A subspace L in X is said to be 
rational if L PI A is a lattice in L. We fix a Lebesgue measure on a rational subspace (L, L n A) as 
above. An affine subspace A is said to be rational if A is a parallel translation of a rational subspace. 
(Note that a rational affine subspace A is allowed to be a translation of a rational subspace by a 
point which is not rational.) For a rational affine subspace A, we fix a Lebesgue measure on A 
which is a translation of the fixed Lebesgue measure on the rational subspace parallel to A. Any 
integration on a subset in a rational affine subspace is performed by using the Lebesgue measure 
normalized in this way. For each vector u € X, let V^j denote the derivative in the direction u. 

For each non-empty subset S in X, let L{S) be the subspace spanned by the vectors y — x 
with x,y & S, which is parallel to the affine hull, denoted by (S), of S. If 5 C Xq, then L{S) is 
a rational subspace in X. Let L be a rational subspace in a rational space (X, A). The natural 
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projection from X onto X/L is denoted by vr^ : X ^ X/L. If L is a subspace in X, let L-*- C X* 
denote the annihilator of L. The quotient space X/L of X by a rational subspace L is again a 
rational space with the lattice 7r/,(A). 

An inner product Q on a rational space A) is said to be rational if Q{x,y) E Q for each 
x,?/ € Xq. Let Q be a rational inner product on (X, A). The rational inner product on {X*,A*) 
induced by the inner product Q on X is also denoted by Q. Let L be a subspace in X. The 
orthogonal complement of L in X is denoted by L-^Q . Note that we have a natural identification 
(X/L)* = L^. The orthogonal projection from X* onto {X/L)* = L"*- is denoted by pL : X* — > 
(X/L)*. When L is rational, the rational space X/L is equipped with the rational inner product 
obtained by identifying X/L with Note that, with this identification, the lattice 'Kl{K) of 

X/L is identified with the orthogonal projection pl(A) of A, where the orthogonal projection from 
X onto -L^<3 is also denoted hy pL : X ^ L^Q , which is different from the lattice L^'^ n A in L^<3. 

A subset P in a rational space (X, A) is called a rational polyhedron if P is an intersection of a 
finite number of half spaces each of which is bounded by a rational affine hyperplane. Let P be a 
rational polyhedron. Then the set of faces of P is denoted by J-{P), and, for non-negative integer 
k, the set of faces of P of dimension k is denoted by F{P)k- We set V{P) = T{P)o, the set of 
vertices of P. A face of codimension one is called a facet. For each / G J^iP), we set vrj = T^L{f)-: 
the natural projection from X onto X/L{f). When, a rational inner product on X is fixed, we set 
Pf = PL{f)j the orthogonal projection from X* onto (X/L(/))*. A rational polyhedron C in X is 
called a rational cone if C is a cone generated by a finite number of elements in A. Note that a 
rational cone C might contain straight lines. The largest subspace contained in the rational cone 
C is C n (— C), which is a rational subspace in X. If C fl (— C) = {0}, then the rational cone C is 
said to be pointed. If a rational cone C is generated by a subset of an integral basis of A, then C is 
said to be unimodular. A subset a of X is called a rational affine cone if a is of the form a = s + C 
where s € Xq and C is a rational cone. If C is pointed, then o is also said to be pointed. 

1.2 The Berline-Vergne operators 

In this subsection, we recall the construction of operators given in |BeVj . Let (X, A) be a rational 
space with a rational inner product Q. For each rational polyhedron P in X, we set 



if the sum and the integral converge absolutely, where ^ G X*. These functions are defined as 
meromorphic functions on X*. Let / be a face of a rational polyhedron P in X. Let Cp{f) be the 
cone generated by the vectors of the form y — x with y P, x € f . This is actually a rational cone 
in X with Cp{f) n (-Cp(/)) = L{f). Then, the pointed affine cone t(P, /) := 7r/(( / ) + Cp{f)) in 
X/L{f) is called the transverse cone of P along /. 

For any rational quotient W = X/L of X by a rational subspace L, let C{W) denote the set 
of all rational affine cones in W . Let T-LiW*) denote the ring of analytic functions with rational 
Taylor coefficients defined in a neighborhood of in W* with respect to an (and hence all) integral 
basis of the dual lattice of the lattice ttl{A) vhW = X/L. 

Then, it is shown in Theorem 20 in |BeV| that there is a unique family of maps ^jlw-, indexed 
by rational quotient spaces W of X, from C{W) to T-LiW*) such that the following conditions hold: 




7GPnA 



(1.1) 



(1) If = {0}, then /ivF({0}) = 1. 



(2) If the affine cone a € CiW) contains a straight line, then ^w{^) = 0. 
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(3) For any a E C{W), one has 



5(a)(0= Yl t^w/L(F){i{a,F)miiFm, ^^W*. (1.2) 

Moreover, one of main theorems in |BeVj is that, for each rational polyhedron PmW = X/L, one 
has 

s{pm= E f^x/LifMPjmwm, ^^w*. (i.s) 

(See Theorem 21 in |BeV] .) Note that the functions Hx/L(f) iii (jl-Sp (and also in (jl.2p ) is the lift to 
W* of functions defined on {W/L{f))* through the orthogonal projection pf : W* {W/L{f))*. 
Let be a pointed rational affine cone in the rational quotient X/L of X. For any non- negative 
integer k, let /i^y^(o) denote the homogeneous polynomial of degree k on {X/L)* which is the 
homogeneous part of the Taylor expansion of the analytic function ^Jix/hi^) ii^ar G {X/L)* . We 
set lJi\{a) = p*^fi^^j^{a), which is a homogeneous polynomial of degree k on X*. 

Definition 1.1 Let (^, A) be a rational space with a rational inner product Q. For any rational 
polyhedron P in X , any face f of P and any non-negative integer n such that n — dim{P)-\-di'ni{f) > 
0, we define the homogeneous differential operator D^{P; f) on X with rational constant coefficients 
of order n — dim(P) + dim(/), which involves derivatives only in directions perpendicular to the 
subspace L{f), as the differential operator whose symbol is given by /U^ '^^'^iP)+^^"^if) ^^(^p^ j^'^ _ 
Pj^^^^™|^^'^'^™^''^^(t(P, /)). We call the operators D^(P;f) the Berline-Vergne operators. 

We note that, when C is a pointed rational cone in X and F is a face of C, then t(C, F) = Trp{C), 
and hence we have D^{C;F) = {-7Tf{C);0). Let P be a lattice polytope in X, which means 
that each vertex is an element in A, and let / G J^{P)- Then, we have t(P, /) = nf{v) + T^f{Cp{f)) 
where v G fOA. Since the function Hx/L{f) is invariant under translation by elements in the lattice 
(Theorem 21 in [BeV] ). we have D^{P;f) = P»^(7r/(Cp(/)); 0). 



1.3 Heuristic arguments 

In this subsection, we give a heuristic argument to find the formula ()0.2p by using the result ()1.3p 
in |BeV| . Let (^, A) be a rational space. Let P be a lattice polytope in X. For simplicity, assume 
that m := dim(P) = dim(X). For each / G F{P), we set fj,{P, f) := p*ftJ'X/L{f){KF^ f)) which is 
a meromorphic function on X* analytic in a neighborhood of the origin. Now let us compute the 
Riemann sum Rn{P; ip) by using (jl.Sp . Let (fhe a smooth function on P. Since P is compact, one 
may assume that (p G Cq°{X). Normalize the Lebesgue measure d^ on X* so that it satisfy the 
Fourier inversion formula 

V9(x) = (27r)— / e^<?'-)<^(Ode, m = [ e'^^^'^^ cp{x). 

JX* JX 

Inserting the above for x = ■y/N with 7 G A^P H A into the definition of Rn{P', ^) and using the 
formula ()1.3p . we have 

Rn{P;^) = ^^^^ J]^^/i(A^P,A^/)(iC/A^)/(A^/)(ie/A^)^(Orf^, 



7 



But, since P is a lattice polytope, we have ij-{NP, Nf) = iJ,{P, f) (see |BeVj . Remark 29). Changing 
the variable x^x/N, we have I{Nf){i^/N) = iV<^™(^)/(/)(i^). Thus we have 

Formally, substituting the Taylor expansion 

k>0 

into the above formula, we could have 

ii^(P;^)"~"5]iV-" Yl /l)^(P;/)^, (1.4) 

where D^{P;f) is defined in Definition 11.11 However, the above computation is formal because 
we do not know much about global properties of the functions n{P, /). Even if we could prove the 
formula (II. 4p along with the method explained above, we do not know much about homogeneous 
parts of its Taylor expansion. One of our purposes in this paper is to give an effective formula for 
the operator D^(P;f) given in Definition 11.11 at least for Delzant lattice polytopes, by a method 
different from the above strategy. 

1.4 A uniqueness property 

In this subsection, we discuss a uniqueness property of an expression of each term of the asymptotic 
expansion of Rn{C;(p) for unimodular cones C. Let C be a unimodular cone in a rational space 
(X, A) with a rational inner product Q. Then, note that, for each face F of C, we have t{C,F) = 
ttf{C). Note also that, we give a rational inner product in each rational quotient space X/L by 
identifying X/L with L^Q . 

Theorem 1.2 Suppose that, for any rational space (-'^, A) with a rational inner product Q, any 
rational subspace L of X, any unimodular cone C in X/L and any non-negative integer n such that 
n > dim(C), there exists a homogeneous differential operator V^{C) on X of order n — dim(C) 
with symbol v^{C) such that 

(1) IfC C X/L, then u^{C) = plvn'^{C) where pL : X ^ L^Q ^ {X/L)* denote the orthogonal 
projection. 

(2) IfCcX with dim(C) < dim(A:), then v^{C) = h*ciyn''^\c), where he : X* L{C)* is 
the transpose of the inclusion lc ■ L{C) ^ X . 

(3) When dim(X) = 0, we have ^^({O}) = 1, P^f ({0}) = (n > 1). When dim(X) = 1 and 
C = M+n with a generator u of K, we have (C) = — ^V"""*^ {ri> 1). 

(4) For any unimodular cone C <Z X, any F € J~iC), any n £ Z+ with dim(F) > dim(C) — n 
and any Schwartz function ip € S{X) on X , the following holds: 

Rr,iC;ip)r.Y,N-^ E j V^{^f{C))p (iV^oo). (1.5) 

n>0 FeT{C);d\m(F)>dim{C)-n^ 
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Then, we have 

'^^{C) = f,-^-'-'-^''\c) (1.6) 
for any such X , C and n satisfying n — dim(C) > 0. 

Proof. First, we note that, the symbols of the Berhne-Vergne operators satisfy the assumption 
(2) in the statement (Proposition 13 in |BeV| ) . 

We prove the assertion by induction on the dimension of X. Consider the case where dimX = 1. 
Take a generator u of the lattice A and identify u with 1 in Z. Let C = M+n. Then, as is computed 
in |BeV| . we have 

1 1 °° t 

^ ' n=l 

We also have /X{o}({0}) = 1. From this, we have /i^"HC)(0 = -^(C,"")""^ (n > 1), fJ-x{{0}) = 1, 
/i^({0}) =0 (n > 1). By the assumption (3), this shows the assertion when dim(X) = 1. 

Next, assume that, for any rational space (^, A) with dim(X) < m — 1, any unimodular cone 
C in a rational quotient X/L and any non- negative integer n such that n > dim(C), the equation 
fll.6p holds. Let X be an m-dimensional rational space, and let C C X be a unimodular cone. 
If dim(C) < m, then by the assumption (2) and the induction hypothesis, we have (jl.6p . Thus, 
we assume that dim(C) = m. Let F S J^{C). If dim(F) > 0, then, by the assumption (1), we 
have {7Tf{C)) = p*pVn^^^^\'Kp{C)). Since d\m{X / L{F)) < m — 1 and irpiC) is a unimodular 
cone in X/L{F), we can use the induction hypothesis, and hence the latter function coincides with 
P*Fl^x~/TtFt'"^^\''F{C)) = /i7™+'^''"(^)(7r^(C)). To prove {C) = ^^^-"(C) for n > m, take 
^ G X* such that {S,,x) < for each x G C. Then, for any > 0, we have 

Note that there is a G S{X) such that ^{x) = e^{x) for x £ C. Thus, by the assumption (4), we 
have 

n>0 FgJ'(C) ;dim(F)>m-n 

as A^ ^ oo. By (fO]) and the identity I{F){i/N) = A^^™(^)/(F)(0, we have 

s{cmN) = N" pT"''"''^^^\Mcmnm) 

n>0 FeT{C);dim{F)>m-n 

for every sufficiently large A^. Let n > m. By using the induction hypothesis, the coefficient of 
j^m-n -j^ ^YiQ above can be written as 

F&T{C) ; 0^dim(F)>m-n 

Equating ^ with the coefficient of A^™"" in ^ shows ^^"""(C) = v^{C). □ 



2 Szasz functions and their asymptotic behavior 

In this section, we define Szasz functions over unimodular cones and investigate their asymptotic 
behavior. First of all, let us compute in one dimension, which illustrate the general case. 
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2.1 Computation in one dimension 

The Szasz function associated with a function (/? on M, originahy introduced and discussed in [S], 
is defined by 

oo ^ 

SN{^){x) = Y,^k{Nx)ip{k/N), 4(x) = |^e-^ x€R. (2.1) 

k=0 

Szasz introduced the function Sj\f{(p) as an analogue of the Bernstein polynomial 
Bn{^){x) = Y,m%(x)ip{k/N), m%{x) = (^^ Ja:''(l " 

Indeed, these two functions are related through Poisson's law of rare events 

lim m%{x/N) = 4(x). 

For us, an important property of the Szasz function S^if) is the following: 

/ SNiv){x)dx = -y2^{k/N) =: i?jv([0,+oo);<^) 
for any ip £ 5(R). We put 

^ oo 

i?iv((-oo,0];(^) := -Y^^{-k/N). 

k=0 

Then, once we obtain the asymptotic expansion of S]sf{(p) as — )• oo with a suitable reminder 
estimate, then integrating it on [0,oo) will give the asymptotic expansion of i?7v([0, +oo); 93). But 
then we have the formula 

Rn{[0, 1]; vp) = Rn(.[0, +00); if) + i?iv((-oo, 0];Tiip) - Rn(R; (2.2) 

where we set Tiip{x) = 99(1 + x). In this formula, note that we have Rj\[{l&;ip) = J^ip{x)dx + 
0{N~^) (see [ GSj or see Lemma [312]) . We also have i?Ar((— 00, 0]; Ti(/?) = i?7v([0, +00); V'), where we 
set tp{x) = (p{l — x), and hence the asymptotics of R]\f{[0, +00); (p) will give the classical asymptotic 
Euler-Maclaurin formula (10. 5p . Thus, to obtain (10. 5p . it is enough to consider Rn {[0, +00); ip). In 
one dimension, we can consider a bit more general situation. We choose a positive integer g > 1 
and a q^^ root of unity u. We consider the twisted Riemann sum 

^ 00 

fc=0 

where (p € Cq°(M). The twisted Riemann sum R%{ip) is discussed in |GS| and the asymptotic 
formula 

RU^) - jVn (2-4) 

n>l 

was obtained, where the coefficients 6^ is defined by the Taylor expansion of the function 

r.{s):= ' ^^ =Y^b-s\ ^i=T^- (2-5) 
1 — we ^ ^-^ 1 — bj 

n>l 
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The formula (|2.4p is used in [GSj to obtain asymptotic Euler-Maclaurin formula for simple poly- 
topes. Now, to obtain the asymptotic expansion of the twisted Riemann sum along with 
our strategy, we use the twisted version of the Szasz function, which is defined by 

oo 

S%{v)ix) =Y,^''^k{NxMk/N). (2.6) 

A;=0 

Prom the definition, we have 

/ S%{ip){x)dx = R%{^). (2.7) 
Jo 

To state a result on asymptotic expansion of the twisted Szasz function S'^{ip), we need to prepare 
some properties of the Stirling numbers of the second kind and related polynomials. 

The Stirling numbers of the second kind, denoted by S{n, k) where n, k are integers satisfying 
< A: < n, are defined by the following recursion formula: 

5(0,0) = 1, 5(n,0) = 0, 5(n,n) = l (n > 1), 

S{n + l,k) = kS{n, k) + 5(n, A: - 1) (1 < A; < n). ^ ' ' 

For example, we have S{n, 1) = 1 (n > 1) and S{n,n — 1) = (2) (n > 2). For convenience, we set 
S{n, k) = for < n < A;. For any integer n, k with < A; < n, we define the polynomial p{n, k; z) 
in z G C of degree k by 



p(n, k- z) := (fj {-^YS{n -t,k- t)z^-'. 



(2.9) 



Some of p{n, k; z) are computed as follows. 

p(0,0;z) = l, p(n,0;z)=0, p(n, n; z) = (z - 1)" (n > 1) 



n 



(2.10) 



p(n,l;z) = z, p{n,n-l]z)=\^;^y{z-lY (n > 2). 
Lemma 2.1 (1) For any non-negative integer n, we have 

n °° i,n 

e^Y.S{n,k)z' = Y.-^zK 

k=0 k=Q 

(2) The polynomials p{n, k; z) satisfy the following recursion formula: 

p{n + 1, k\ z) = [z — l)p{n, k — 1; z) + kp{n, k; z) + np{n — l,k — 1; z), 1 < k < n. 

(3) For [n/2] + 1 < k < n, the polynomial p{n, k; z) is divisible by {z — 1)2'^-"-, particular, we 
have p{n, A;; 1) = for [n/2] + 1 < k < n. 

Proof. (1) is proved easily by using induction on n and the recurrence formula for the Stirling 
numbers S{n,k) of the second kind. To prove (2), let 1 < A: < n. By using the relation ("+^) = 
(t) + it-i) for 1 < t < n, we have 



p{n + l,k;z) = Y^ (^^ (-l)*'S'(n + l-t,k-t)^ 



z''"* -p{n,k- l;z). 
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Denote S the sum above. Then, by the recursion formula (|2.8|) . we have 

k . ^ fc-l 



^ = i2(1) - *)^(^ - ^ - *)^'"* + E (?) -t,k-l- t)z''-' 

= kp{n,k;z) - (j^ y\ {-if S{n - t,k - t)z^~^ + zp{n,k - l;z). 
t=i ~ ^ 

Minus the sum in the middle of the above equals np[n — 1, A: — 1; z), and hence (2) is proved. 

Let us prove (3). Since the statement is obvious from (j2.10p for n = 1,2, we assume that, for 
some n > 2, p(m, k\ z) is divisible by (z — fg^. g^ch 1 < m < n and [m/2] + 1 < k < m, 

and use the induction on n. So, we take / with [(n + l)/2] + 1 < / < n + 1. If / = n + 1, 
p{n + 1, n + 1; z) = [z — 1)""^^ and hence (3) is clear. Thus, we assume that [(n + l)/2] + 1 < / < n. 
By the induction hypothesis, p{n, I; z) is divisible by (z-l)^'-". We have [(n-l)/2] + l = [(n+l)/2] 
and hence, by induction hypothesis, p{n — 1,1 — 1; z) is divisible by (z — i)2'-n-i^ jf =1 — 1, 
then n is even and 2/ — n — 1 = 1, and hence, by the recurrence relation (2), p{n + l, I; z) is divisible 
by (z — 1). Otherwise, we have [n/2] + 1 < l — l, and hence p{n, l — l;z) is divisible by (z — i)2/-n-2_ 
Then, again by (2), p{n + 1, /; z) is divisible by {z — i)2/-n-i_ \^ 

Now, we can state the asymptotic expansion of the twisted Szasz functions S'^^cp) by using the 
polynomials p{n, k; z) as follows. 

Proposition 2.2 Let Lp E Let u be a g*^ root of unity. Then, for any positive integer n and 

positive number K such that n < K < 2n, there exists a constant CK,n > such that we have 



2ra-l 



SUv){x) = ^^^N-f^J^iNx) + S^^^j,{x), X > 0, (2.11) 

where the function iy{x) satisfies the following estimate: 

\S^nM'')\ ^ CK,nN~''{l + x)"-^, X > 0, iV > 0. (2.12) 
The function J^{x) is given by 

r^{x) = e-^^-^^^j2p{^^,k■uJ)x\ (2.13) 

fc=0 

When uj = 1, the function Jj^{x) is a polynomial in x of degree at most [fJ-/2]. 
Proof Let x > 0. Substituting the Taylor expansion 



^{k/N) = ^H^ik/N - xY + R2n{k/N,x), 

^-^ u! 2n — 1 ! 

0</i<2n-l ^ ^ ' 

R2n{k/N,x) = C {l-tf''~^ip^'^"'\x + t{k/N -x))dt. 



we have 
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where and S2n,N{x) are given by 

oo 

J-(x)=^a;'=4(x)(fe-xr, 

k=0 

_ oo 

S^nM^) = - i)!Ar2n Yl^%iNx){k - Nxf^R2n{k/N,x). 

By using Lemma [2 .11 (1) and the definition (|2.ip of the function ikix), it is easy to show the formula 
(|2.13p for J^{x). We set 82^^ ^{x) = ^2n-i)\N'^" S2n,N{x)- Take K as in the statement, and choose 
C > so that [v7(^"Hy)| < C(l + |y|)--^ for any y £R. Then, we have \x + t{k/N-x)\ > {l-t)x for 
any t G [0, 1], x > 0, A; > 0, and hence \R2nik/N,x)\ < CK,nX~^, x > I, k > 0. When < x < 1, 
\R2n(k/N,x)\ is bounded uniformly in A''. Thus, we have |5'2n,Af(a;)| < Cx~^ j\^{Nx) for x > 1. 
When < X < 1, we have |52„,Ar(x)| < CN-'^'^Jl^{Nx). But, by LemmaEH (3) and the formula 
()2.13p . J2n{x) is a polynomial in x of degree at most n. Therefore, we obtain (|2.12p . □ 

In general, for any r G C with Re (r) > and any n > 0, we have 

Taking > in Proposition 12.21 so that n + \ < K <2n and integrating ()2.1ip . we conclude the 
following. 

Proposition 2.3 When lo ^ 1 is the g*^ root of unity, we have 

c„ = ^^ ■ (2-14) 

n>l o=Ofc=0 ^ J \ / 



When oj = 1, we have 

RN{[0,c^y,v)^ / vix)dx + ^c j c„ = ^i^Zi^(-l)-«+Va,a-n), (2.15) 

•^^ n>l a=n ^' 



where we set 



p(n,fc) :=p(n,A;;l) = ^ (''^V-l)*5(n-t,A;-t), < k < n. (2.16) 

Note that a direct computation and the well-known formula for the relation among the Bernoulli 
numbers, Catalan numbers and the Stirling numbers ( [GKPj ) shows 

(n+l)^2n\-i " (-l)Y 2n \ 



n! \ n , 
which shows that, for w = 1, we have 



RN{[0,oo);ip) ^ ^(x)dx-y^^ip^"-'\0)N~^, (2.18) 

from which we have ()0.5p . For 7^ 1, we compare each term of the asymptotics (j2.4p . (j2.14p to get 

6^ = (-l)«-ic^ = V y:(-l)^+^ (n-fc-1)! p(«,a-fc;^) ^ 
n V ; n Z^Z^v ; Q;!(n-Q;-1! w - 1 ""'^ ^ ^ 

a=0/c=0 ^ J \ J 
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2.2 Definition of Szasz functions 

Let C be a unimodular cone in X. Since the Riemann sum Rn{C; ip) depends only on the restriction 
of if to L{C), replacing (X, A) by {L{C),L{C) n A) if necessary, we assume, for a moment, that 
dim(C) = dim(X). Then, C is written in the form 

e&E 

where E is an integral basis of A, and R+ denotes the set of non-negative real numbers. For abstract 
two sets S and T, let S'^ be the set of all functions from T to S. The whole space X is identified 
with M^. Since E is an integral basis, A is identified with Z^. Note that, C and CnA are identified 
with and Z;^, respectively, where denotes the set of non-negative integers. For any a € 
and X G AT, we set 

a! = JJa(e)!, = JJ x(e)"(<=), 

e&E eeE 

where x{e) is the value oi x at e G E when we identify X = M^. For each 7 G Z;^, we define the 
function on X by 



7! 



.e-EeeB^(e). (2.20) 
Then, the function is non-negative, integrable on C and satisfies 

f i^{x)dx = l ^iy{x) = l (xeX). (2.21) 



Definition 2.4 We define the Szasz measure S{x) = dSx on X = M^, parametrized by x ^ X , by 

S{x) = dSx := t^{x)5r^ 



^1- 



By the second property of (j2.2ip . the measure dS^ is a probability measure on C. For each G 
the A^-th dilated convolution powers, denoted by dS^ , of dSx is given by 

dS^ := {DyMS{x) * • • • * S{x)) = ^iiNx)6^/N, 
where Di/^ : X ^ X is the dilation Di/^{x) = x/N, x G X. 

Definition 2.5 We define the Szasz function SNif) associated to a function ip on X, by 



Sn^x) := I ^iz)dS^{z) = ^ £^iNxM-f/N) (2.22) 



if the sum in the right hand side converges absolutely. 
By ()2.2ip . the Szasz function Sn{'p) satisfies that 



^^(^;^)^=]^dic) E V>h/N) = I^SN{^){x)dx (2.23) 



7GCnA 

if the sum converges absolutely. 
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2.3 Asymptotics of Szasz functions 

For each n^v & with fj, > we define 

PEif^,i^) = llpiKe),'yie)), (2.24) 

eeE 

where p{n,k) is an integer defined by (12.16p . For each fi € Z;^, we set = Ilee-B ^^'^ 
\fi\ = X]eg£;/^(^)- Then, a relevant asymptotic formula for the Szasz function Sjyi^) is given as 
follows. 

Proposition 2.6 For each positive integer r and positive number K with r < K < 2r, there exists 
a positive constant Cr.K such that we have 



Sn{^){x)= -j^MNx) + S2rM^), (2.25) 

^teZf ; |/^|<2r-l ^' 



where the function S2r,Nix) satisfies the following estimate; 

|^2r,7v(^)| < Cr,KN-''il + \x\Y'^ , X G C, (2.26) 

where the norm \x\ ofx ^ X is defined by [xp = J2e€E ^i^)"^! '^'^^ the function J^{x) is a polynomial 
in x of degree at most [\fi\/2] given by 

Mx)= Yl Pe{^i.v)x^. (2.27) 
where [/i/2] G Zf is defined by [^i/2]{e) = [^i{e)/2]. 

Proof. The proof is the same as that for Proposition 12.21 Inserting the Taylor expansion 

^ ZM^(,_^)M+ Y ^R2rA^,x)iz-xr, 

; |^t|<2r-l lAt|=2r- 

-R2r,/.(-z, x) = 2r f (1 - tf^-^Vifix + t{z - x)) dt 
Jo 

with z = j/N into the definition (|2.22p of the Szasz function Sj\f{(p), we have 

where the functions J^(x), S2t,^i{x) are given by 

J^{x)= Y iyiNx)i^-xr, 

S2rAx)= Yl 5] ^(iVx)i?2.,;.(7/iV,x)(7-iVx)^ 
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The formula (|2.27|) is easily obtained by the relation 

which follows from Lemma |2.H (1), where Se{i^iO) is given by 

SE{u,a) = \{S{v{e),a{e)). (2.28) 

Next, we estimate the term S2r,N{x)- Note that x and ^/N are in C . Thus, we have \x + t{'y/N — 
x)\ > (1 — t)\x\ for each < t < 1. We choose a positive constant Cr^K such that \V^ip{x)\ < 
Cr,K{^ + l^^l)^^ for each /i G with = 2r. Hence, if |x[ > 1 and = 2r, we have 

\V^ip{x + t{-f/N - x))| < Cr,K{l + \x + t{-f/N - x)\)-^ < Cr,K{l - t)-^\x\-^ (2.29) 

Thus, for |x| > 1, we have \R2r,fii'y/N,x)\ < Cr,K\x\^^ < Cr,K{'^ + \x\)~^ , where Cr,K is a constant. 
Therefore, we obtain 

\S2rM^)\ < ^(1 + Ixi)-^ Yl ^(^^) E - ^ §#(1 + l^i)""" E -^^mIa^^)- 

7GZf lMl=2r- ^' |At|=r 

As is mentioned above, the function J2^(x) with = r is a polynomial in x of degree at most r. 
Thus, we have | J/i(x)| < C^|a;|^ where does not depend on x. Therefore, we obtain the estimate 
(12.26j) . When < 1, we estimate R2r,^i{l/N,x) as \R2r,ij.{l ,x)\ < Cr,K, and hence S2r,N{x) is 
bounded by Cr,KN~^ < Cr,KN~^{l + Ixl)*""^, which completes the proof. □ 



3 Asymptotic Euler-Maclaurin formula over unimodular cones 

In this section, we deduce asymptotic Euler-Maclaurin formula of the Riemann sum over unimodular 
cones in a rational space {X,A). At first, we deduce it by using Proposition 12.61 The result 
coincide a well-known result due to Guillemin-Sternberg |GSj . We don't need to use a rational 
inner product on X so far. Then, we renormalize each term of the expansion using an integration 
by parts procedure to find explicit form of Berline-Vergne operators. This step involves a rational 
inner product. 

3.1 An Euler-Maclaurin formula for unimodular cones 

As before, let C be a unimodular cone in X with dim(C) = dim{X) and let E be the integral basis 
of A generating C. For each I C E, let |/| be the number of elements in /. For such /, we regard 
Z;?^ as a subset of consisting of a G with the property that a{e) = for each e G E \ I. 
Clearly we have C Z'^ if I C J. For any e G E, we define Ae G by Ae(e) = 1, Xe{v) = 0, 
V G E \ {e}. Then, we obviously have a = Ylie&E 

(e)Ae for each q G . For / / C E, we set 
^>o = {" G ^+ ; oi{e) 7^ 0, e G /}. For / = 0, we set Z^q = {0}. Each I C E corresponds to a face 
C{I) of C defined by 

C{I) := Y ^+^' ^(^) ■= {0}, (3.1) 

eeE\I 
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and for each face F of C, there is a unique Ip C E such that F = C{Ip). Thus, we identify subsets 
in E and faces of C. Note that F C G if and only if Iq C If- For each /i, € with u < fi and 
/ / C we set 

:=nP(Me),Ke)). (3.2) 

If ^, G we clearly have pj{^,v) = pi{^,v) for each J with I C J because p(0, 0) = 1. For 
each u G Z;^, we set 

P/H= E (3.3) 

Then, we have pj{i') = pi{i^) if v ^ Tllj^ and I d J . Note that pi{v) = Y\e(^iPi^i^))^ where we have 
p{n) = {-l^-^Cn = (-l)"6„/n! as in (I2T5I1 . (fZTTD . 

For each non- negative integer n and a subset I of E with |/| < n, we define a homogeneous 
differential operator Ln{C;I) of order n — \I\ on X with constant coefficients by 

L„(C;/) = (-!)" E P/(^)V"-^('\ e(/) = EAe (n > 1), (3.4) 

and Lo(C; 0) = 1. When n > 1 we set 1^(0; /) = for |/| > n or / = 0. 
Proposition 3.1 For each if E S{X), we have 

RNiC;^)^^N-" E (-1)'"/ LniC;!)^. (3.5) 



n>0 /C-E;|/|<n 



Proof. We use Proposition 12.61 We take r € N and > so that r + dim(X) < K < 2r. By the 
estimate (j2.26p . one can integrate the asymptotic expansion ()2.25p over C. Then, by ()2.23p and 
([2:25]) . we have 

Rn{C;^)= E ^N-\^^-'\pE{^i,u) j x'^V^'^ + 0{N-^). 

fi,u;\fi\<2r-l,u<[fi/2] ^' 

Integrating by parts, we have / x'^V^ip = (— l)''^'!^! / \7^~'^ip, and hence, substituting this into 

Jc Jc 
the formula for R]\f{C;ip) above, we obtain 



r-l 



(3.6) 

Lk{c) = i-i)' E PEi^)'^''= E (-l)l^l+'=^^^^PE(^,^-^)v^ 



u£Zf,\u\=k u,ii,\u\=k,u<fj.<2u 



To integrate by parts further in the right hand side, note that we have Z^ = IJj^^jZ^q, which is 
a disjoint union. For u G Z{,q, we have 



Jc Jc{i) 
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From this, we obtain 

/ L,iC)^= Yl (-1)'" / L,{C-I)ip, 

''^ I<ZE;\I\<k 

which shows the assertion. □ 

Next, we consider cones containing straight hnes. Let E be an integral basis of A, and let 
I C E. Consider the cone C in X of the form 

C = ^M+e + L, (3.7) 

where L is a subspace in X spanned by vectors e (z E\I. If I = E, then L = {0} and in this case 
C is a unimodular cone discussed above. When / = 0, we set C = X. 

Lemma 3.2 Let E be an integral basis of A, and let C be a cone of the form ()3.7p with I C E. 
Then, for each E C^{X), we have 

Rn{C; v9) = Rn{7Tl{C); (ttl)*^) + 0(iV-°°), (3.8) 

where {'Kl)*'-P is a compactly supported smooth function on X/L defined by 



/ (p, X e X/L. 



Proof. For simplicity, we write vr = ttl : X ^ X/L for the natural projection. Take ip € C^{X). 
For any G X, we set TyLp{x) := ip{x + v). Let M be the subspace spanned by / so that X = M®L. 
We identify L with M^\^ and M with in a natural way. Then, we can choose v G Z^^^ so that 
supp(T„(/?) C M + R^y. Clearly we have Rn{C;(p) = RnO^^'tT^^), where we note that R;^ is a 
unimodular cone in X. Therefore, by (I3.6p . we have 



Rn{C; v?) ~ V N-^ [ Ln{Rf)nip, 

where the differential operator L„(R^) is given in ()3.6p . Note that 7r(C) is a unimodular cone in 
X/L with respect to the lattice vr(A) generated by the integral basis '7r(/) of vr(A). Since v S Z^'^^, 
we have 7r*T„(^ = -K^cp. Therefore, according to Proposition 13.11 we only need to show that 

Ln{Rf)^pdx= / Lni-KiC))-K^7pdx (3.9) 

for any ^ G Cg^i^) with supp (ip) C M + R^^. If i/ € has some e e E\L such that z^(e) > 1, 

then, since supp (ip) n R^^^*^-*^ = 0, we have / V^tp = and hence 

Jmf 

f ^Ln{Rf)iPdx= f ^Ln^P, Ln = {-ir Yl Pii^)^"- 

If we denote = Heg/ ^!^(e) each v £ L, then, by the definition of the function Tr^ip on X/L, 
we have V^vr^.V' = Tr^V^tp for each € Z;'^. Since supp (ip) C M + R^^^, we obtain, for u S 

7r(C) J7r(C) ^Kf 
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From this and the definition of L„(7r(C)), we obtain (|3.9|) . 



□ 



Remark: As is mentioned in Introduction, Proposition 13.11 is deduced directly from Theorem 3.2 
in [GSj. Lemma [32] is also obtained in [GS] . 

3.2 Integration by parts 

In Proposition l3.1l and Lemma r3.2l in the previous subsection, we have derived an asymptotic formula 
for the Riemann sums over unimodular cones and their variants. In each term in these asymptotic 
formulas, integration over faces of homogeneous differential operators L„(C;/) defined in ()3.6p 
appears. The differential operators L„(C;I) involve derivatives only in directions transversal to 
the face C{I). However, these derivatives are not 'perpendicular' to the face C{I), and hence we 
can perform further integration by parts. If one performs integration by parts in (j3.5p . then one 
will find the differential operators which involves derivatives only in directions perpendicular to 
faces. However, we need to perform this procedure systematically to define the operators all at 
once. This step is one of the main points which makes the final formula complicated. 

In the rest of this paper, we fix a rational inner product Q on the rational space {X,A). Let E 
be an integral basis of A. For each I C E, we set 

X{I) = Me ^ M^\^ X{E) = {0}. (3.10) 

e&E\I 

Note that X = X{$), and if I C J, then X{J) C X{I) and hence X{I)^Q C X{J)^Q. As before, 
for each I C E, we define the unimodular cone C{I) in X by (13. ip . For each a € Z;^, we set 

Proposition 3.3 There exists a family 

{L{E;I,J;a); $ I C J C E,a e Z{,\J\ < |a| + |/|} 

of homogeneous differential operators L{E; I, J; q) of order \a\ — \J\-\-\I\ on X with rational constant 
coefficients which involves derivatives only in directions perpendicular to the rational suhspace X{J) 
such that for each {I, a) with / / C -B, a € 1^^^, we have 

/ V>= Yl (-l)l^l-l'l / L{E;I,J;a)cp (3.11) 

'^^^^ J;/cJ,|J|<|a|+|/| 

for any ip € SiX). Furthermore, fix a and ^ ^ I (Z E with a G 'L^_^. Suppose that a family 
{L{J) ; / C J C -E, I J| < |a| + |/|} of homogeneous differential operators with constant coefficients 
of order \a\ — \ J\ + \I\ which involves derivatives only in directions perpendicular to X(J) satisfy 
the equation (j3.1ip for any ip € S{X). Then, we have L{J) = L{E;I,J;a). 

Note that a differential operator on X is said to have rational coefficients if it has rational 
coefficients with respect to an (and hence all) integral basis of A. We first prove the existence of 
such family of differential operators. 

Proof of the existence in Proposition 13.31 For a given 7^ / C -E and e G /, we decompose e along 
with the orthogonal decomposition X = X{I)-^Q © ^{1)7 which is denoted by 

e = n(/;e)+ ^ c(/;e,f)f, u{T,e) G X{I)^Q n Xq, c(/;e,t;)eQ, e G /. (3.12) 

v(^E\I 
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We set u{E; e) = e for each e G E. We construct the operators L{E; /, J; a) inductively as fohows. 

(0) When \a\ = 0, then \ J\ = \I\ and I C J imphes I = J. In this case, we set 

L{E;IJ;0) = 1. (3.13) 

(1) We take % ^ I d J d E and a € with |a| = 1 and \ J\ < \a\ + |/|. In this case, J = I or 
J = I U {v} with V e E\I and a = Xe with eel. We then define L{E; I, J; Ae) by 



(3.14) 



L{E;I,r,Xe) =Vu(i-e) (when/=J), 
L(£;; 1,1 U {v}; Xe) = c{I; e, v) (when J = / U } with v e E\I). 

Then, by the identity 

/ V,(^ = -/ y,, veE\l, ^eC^iX), 

Jc{i) Jc{IU{v}) 
it is easy to show that the operators L{E; 1,1 U {v}; Xe) satisfy 

/ Ve9?= / L{E;I,I;Xe)ip+ {-!) [ L{E; I, I U {v}; Xe)ip. (3.15) 

JC{I) JC{I) ^g^^^ Jc{IVJ{v}) 

(2) Suppose that, for a positive integer n > 2, we have defined differential operators L{E; I, J; /3) 
satisfying the formula (ISTTTl for any ip G C^(X) for each /, J, /3 with fj) ^ I C J C E, /S e Z{ 
satisfying \/3\ < n - 1, \ J\ < \p\ + \I\. 

(3) For % ^ I d J d E and a € satisfying |q| = n and [ J| < n + [/|, we take e G / such that 
a{e) > 1. Then, we can decompose a as 

a = Ae + /3, /3gZ^, |/3|=n-l. (3.16) 

We define L{E; I, J; a) by the formula 



L{E;I,J;a) 



( L{E;I,I;P)V^(^l,e) (when J = /), 

L{E; I, J; /3)Vuii-e) + EveJ\i ^U; e, / U {z;}, J; /3) 

(when |/| + 1 < I J| < |/| + \a\ - 1), 
[ Ev(,j\i e, v)L{E; I U J; /3) (when | J| = |/| + \a\ 



(3.17) 



A direct computation shows that the differential operators L{E;I,J;a) satisfy (j3.1ip . □ 

Next, we proceed to prove the uniqueness of such family {L[E; I, J; a)}. Fix a and ^ ^ I d E 
such that a E l^i^. If a = 0, then clearly the operator L{E; I, I;0) satisfying (j3.1ip is just the 
constant 1. If / = E, then for any a G Z^, the operator L{E; E, E;a) is uniquely determined as 
L{E; E, E; a) = V". Thus, in the following, we assume a 7^ and % ^ I d E. 

If the homogeneous differential operators {-Z^(«/)}/cJc-B; |j|<|q|+|7| with constant coefficients of 
order \a\ — \J\ + [/[ which involves derivatives only in directions perpendicular to X{J) satisfy the 
equation (j3.1ip for any ip G S{X), then their symbols a{L{J)) must satisfy the equation 

r= E cj{L{jm) W{i,e), e = \{{M^^'\ 

J;/CJ,|J|<|a| + |/| eeJ\7 e&E (3.18) 

a{L{J)){i)=a{L{J)){pM)), i d X\ 
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where the symbol cr{D) of a differential operator D on X (with constant coefficients) is a polynomial 
function on X* characterized by cr(L')(^) = e^^De^, = e^^'^\ x G X, ^ G X* . In ()3.18p . 

Pj denotes the orthogonal projection from X* onto the annihilator X{J)^ of X[J). Therefore, to 
prove the uniqueness in the statement of Proposition I3.3| it is enough to show the uniqueness of the 
family of homogeneous polynomials {(t(-L(J))} satisfying (IS.lSp . First of all, consider the following 
expression. 

a{IJ-am=pi{ir, (3.19) 



a[I,J;a){^) = '—^ — — — ^ , I C J,\J\ = k + \I\ 

J-ICJ,\J\ = \I\+i eGJ\/ 



(3.20) 



where k is an integer satisfying 1 < A: < |q:|. Note that (7o(/;a) = a{I,I]a) = pi{C)°' is a well- 
defined homogeneous polynomial of degree \a\ on X* . Thus, the above equations (I3.19p . (I3.20p 
define rational functions cr(I, J; q), <Ti(I; a) for 7^ / C J C -E, a € Z;^, \ J\ < |a| + |/|, < i < |a|, 
which are homogeneous of degree \a\ — \J\ + |/|, |a|, respectively. Note also that the functions 
(t(/, J; a) satisfy the second line of (I3.18p . 

Lemma 3.4 The functions defined by (j3.19p . (|3.20p are homogeneous polynomials. 

Proof. First of all, let us examine the function a{I, J; a) with I C J, | J| = 1 + |/|. In this case, 
we can write J = I L) {u} with some u € E\I. By (I3.20p . we have 

o-(/,/U {u};a)(0 ^ 



Take ^ G X{I U {n})"*-, which means that (^,e) = for each e ^ E \ I, e u. Thus, there 
exists an r/ G X{I U {li})"*" perpendicular to X{I)-^ with respect to the inner product Q such that 
QiViV) = 1- Note that {r],u) 7^ 0. Let q(^) = Q{S,.,rf}r] denote the orthogonal projection onto the 
one-dimensional subspace Mr/. Then, we have Pi\j{u} = Pi + I1 and hence 

a{I,I U {u};a){^) = 



which is a homogeneous polynomial of degree |q| — 1. Next, to use the induction, suppose that 
the functions a{I,J;a) with / C J, |J| < A; + |/[ (1 < A; < \a\ — 1) are homogeneous poly- 
nomials of degree \a\ — \ J\ + |/|. By (j3.20p . the functions ai{I]a) {i = 0,...,k) are homoge- 
neous polynomials of degree |a|. Take Jq C E such that / C Jo, |Jo| = A; + l + [/[. Set 
/(O = PJoiO" ~ Z]i=i ck)(p^o(0)- Note that the polynomial / is determined on X{.Jq)^. 
So, let ^ G X{Jo)^. Assume that (C,e) = for some e G Jq \ /, which means that ^ G X{K)^ 
with K = Jo\ {u}. Note that I C K C J, \K\ = k + Since ^ G X{K)-^ C X{Jo)^, we have 
PJ(,(^) = ^, and hence by (I3.20p . 

fe-i 

ak{I-am = a{I,K;am J] ( « ) = T - ^ ^^(^^ ")(^)- 

e£K\I i=0 

This shows that /(^) = for G X{K)^. Thus, the homogeneous polynomial /(,^) is divisible by 
the linear function {pjq{C)^^) for each e G Jq\I. Note that the elements in Jq \ / are linearly 
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independent. Therefore, the homogeneous polynomial /(^) is divisible by Y\eeJo\iiPJoiO^^)^ 
hence a{I, Jq; a) is a homogeneous polynomial. This completes the proof. □ 

Proof of the uniqueness in Proposition 13.31 Fix a, I such that ^ a € Z^, % ^ I C. E. Suppose 
that the set of functions {•s(^)}/cJc-E,lJ|<|a|+|/|> where each s(J) is a homogeneous function on 
X* of degree |q| — \ J\ + [/[, satisfy the equation ()3.18p . Let a{I,J;a) denote the homogeneous 
polynomials defined by (I3TT9]) . IKTOh . We need to prove s(J) = a{I,J;a). Let ^ G X{I)^. Then, 
we have (C,e) = for each e £ E\I. Thus, by (|3.18p . we have s{I){^) = a{I,I]a){^) for each 
^ G X{I)-^. Since s(/) satisfies the second equation of (I3.18j) . we have s{I) = a{I,I;a) on X*. 
Next, take Jq C E such that / C Jq, | JqI = 1 + 1-^1 < |a| + l-^l- We write Jq = I U {u}. Take 
^ G X{Jq)^. Then (^,e) = for each e G E \ Jq, and hence the function s{Jo) must satisfy 

r = fTo(/;a)(0 + s{jQm{^,n), ^ G X(Jo)^. 
Since the function s{Jq) satisfies the second line of (I3.18p . we have 

where q is the orthogonal projection onto the one-dimensional subspace of X{Jq)-^ perpendicular 
to X{I)-^. This equation shows s{Jq) = cr(/, Jo;a) for Jq = / U {u}. Now, suppose that for any 
J C E with / C J, I J| < A; + |/|, A; + 1 < |q|, we have s{J) = a{I, J; a). Take Jq <Z E with / C Jq, 
\Jq\ = k A- 1 + \I\ < \a\ + Take ^ G X{Jq)-^. Since (C,e) = for each e £ E \ Jq, the equation 
(l3l8D shows 

C = Sk{0 + s{Jom U "^(^^= E ^T(/,J;a)(C) n (^'^)' (3-21) 

eGJo\-f -fCJCJo, |J|<fc+|/| e£j\I 

for each ^ G X(Jo)-'-. If J C with I C J satisfy J \ Jq / 0, then we have YleeJ\A^'^) ~ ^ 
and hence Sfc(^) = ^,t=o '^«(^' '^)(?)- Since s(Jo) satisfies the second line of (|3.18p . we have s{Jq) = 
a{I,Jo;a). □ 

Remark: One can prove directly that the polynomials ct(/, J; a) defined by ()3.19p . (|3.20p actually 
a solution to (I3.18p . However, its proof is similar to the proof of the existence in Proposition 13. 3( 
and hence we omit it. 

In the above, we have defined the differential operators L{E;I,J;a) for each ^ ^ I <Z J <Z E 
and a G satisfying |J| < |a| + |/|. But we need to work on the quotient space X/L and the 
unimodular cone '7r(C) which is the image of a unimodular cone C under the natural projection 
■K : X ^ X/L where L is a subspace spanned by a subset of E. To state the next lemma, we 
need to fix some notation. Let E be an integral basis of A. For $ ^ K d E, wq set, as before, 
X{K) = 0^g£;^^Mt'. Let t:k '■ X X/X{K) be the natural projection. For each e G E, we set 
e = 7ri<-(e). Then, we have irxiE) = TTxiK) = {e ; e G K}, and the set -KxiK) is an integral basis 
of the lattice TTxi^) in X/X(K). Note that ttk is a bijective map from K onto t^{E). For each 
^ ^ I C K and a G Ij^^, denote vrx(a) G Zi^^^ the Z^^-valued function on iTKil) defined by 

7ri^(a)(e) := a(e), eel. (3.22) 
We note that, for each a G there is a unique a G with the property that 7ri^(a) = a. 
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Lemma 3.5 In the notation as above, we identify X/X{K) with X{K)-^Q to give X/X{K) the 
inner product induced by the inner product Q on X. Then, for each % ^ I d J d K and a € 
with \ J\ < \a\ + [/[, we have 

LiE;I,J;a) = L{TTK{K);TTKiI),TTK{jy,7rKia)) (3.23) 

where the operator L{7rK{K);TTK{I),'^K{J)]'^K{o:)) is regarded as an operator on X by the identi- 
fication X/X{K) ^ X{K)^Q. 

Proof. For simplicity, we set vr = -kk- Let axil, J', a) denote the symbol of the differential 
operator L(7r(Er); 7r(/), 7r( J); 7r(a)) which is a homogeneous polynomial on {X/X{K))*, and note 
that we have the identification {X/X{K))* = X{K)-^ under the transpose *7r : {X/X{K))* X* 
of vr. Then, the symbol of the lift of the differential operator L{i:{K); 7r(/), 7r( J); vr(Q)) is given by 
aK{1 1 J',(y){pK{C))-, C & X* . The symbols of the operators L{E; I, J;a) for J C K, which are as 
above denoted by a [I, J; a), are determined on X{K)^. By (j3.18p . we have 

r= n ^ex(if)^. (3.24) 

J:ICJCK e&J\I 
kl<l"l + l^l 

In Proposition 13.31 we can replace X by X{E \ K) which is identified, as a rational space, with 
X/X{K). With this identification, the symbols ctk{I, J] ot) also satisfy the equation (I3.24p . Noting 
that the equation (j3.24p is nothing but the equation (|3.18p on X{E\K), and using the uniqueness 
in Proposition 13.31 we conclude the assertion. □ 



3.3 Berline-Vergne operators over unimodular cones 

We use the results obtained in the previous subsections and Theorem II. 21 to find an explicit expres- 
sion of Berline-Vergne operators for unimodular cones. 

Definition 3.6 (1) Let C be a unimodular cone in a rational space {X, A) with a rational inner 
product Q. Assume that dim(C) = dim(X). Let E be the integral basis of A generating C . For each 
F € J~iC), we take, as before, a unique Ip C E such that F = C{If). Then, for each F E J~iC) 
and n G Z+ with dim(i^) > dim(C) — n, we define a homogeneous differential operator 'D^{C;F) 
of order n — dim(C) + dim(i<') with rational constant coefficients which involves derivatives only in 
directions perpendicular to the face F by 

P^(C;F):=(-l)«-^^-(^)+<ii-(^) ^ ^ pj^^)L{E-I,lF-u-e{I)), (3.25) 

IC^Ip uGZlg, \u\=n 

and V^{C; C) := 1, V^{C; C) := (n > 1). 

(2) Let C C X be a unimodular cone. For any F G J^{C) and n € Z+ with n — dim(C) + dim(F), 
let V^{C;F) be the differential operator Vn^'~^\C;F) regarded as an operator on X through the 
inclusion ■ L{C) ^ X, where the operator Vn {C;F) is defined as in (1) replacing {X,A) by 
(L(C),L(C)nA). 

For unimodular cones C in X with dim(C) < dim(X), the differential operator T>^[C;F) is 
characterized by the identity t*cV^ {C; F)(p = Vn^^\C; F)l*c(P for ip € C°°{X). Thus, a direct 
computation using Definition 13.61 (|3.1ip . ()3.4p combined with Proposition 13. II (replacing {X,A) by 
{L{C), L[C) n A) if necessary) shows the following. 
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Theorem 3.7 Let C be a unimodular cone in the rational space {X^h) with a rational inner 
product. Then, for any (/? € S{X), we have 

n>0 FeJ"(C),dim(F)>dim(C)-n 

Let C be a unimodular cone in the rational space (X, A), and let F € J'{C). The order of the 
differential operator T>^[C;F) is n — (dim(C) — dim(i^)), and which is equal to the order of the 

differential operator I>„ (7ri;'(C); 0). Moreover, we have the following. 

Lemma 3.8 Let C be a unimodular cone in (X,A), and let F € J~iC). Then, the operator 

T>^{C;F) coincides with the lift of the operator T)n^^^^\'KF{C);Q) on X through the identifica- 
tion X/L{F) ^ L{F)^Q. 

Proof. Let K be the integral basis of L[C) PI A generating C, and let F = C{Lf) with a subset 
If of K. Then, the cone ■7Tf{C) in the rational space {X/L{F),ttf{A)) is a unimodular cone with 
the generator If = {e;e = 7ri?(e), e G If}- Thus, by Definition 13.61 we have 

p5^(^)(^^(C);0) = (_i)"-dim(^HC)) ^ ^ pj{v)L(T^;I,I^;V-e(l)). 

IcT^VGZlg ■ \v\=n 

The subsets I oi K correspond to the subsets / of If by the projection vr^, and the elements V in 
corresponds to the elements v in 1^1^. Therefore, Lemma 13.51 shows 

L(7^; 1,7^; 17 - e(T)) = L{t:f{If)\ 7rir(/), 7r^(/f ); ttfIi^ - e(/))) = L{K- /, If; v - e(I)) 

as an operator on X. From this, the assertion follows. □ 

Example: In one dimension, it is easy to compute the differential operators V^{C;F). Let X be 
a 1-dimensional vector space with the lattice A. Let u G A be a generator and set C = M+u. The 
faces of C are and C itself. Then, E = {u}. By definition, we have V^{C; C) = 1, V^(C; C) = 
(n > 1). By (j3.17p . we have L{E; {u}, {u}; k) = for k € Z+. Thus, by Definition 13.61 we have 

P^(C;0) = (-l)"-Vn)Vri = -^Vr\ n > 1, (3.26) 

n! 

and its symbol is given by —^(^,u)"'~^, and we have Rn(C; f)^ / <^ — / -^V"~^(/?(0). 

Theorem 3.9 For each unimodular cone C in a rational space (X,A), each face F € ^{C) and 
each non-negative integer n such that dim(F) > dim(C) — n, we have 

V^{C;F) = D^{C;F), 

where T>^{C; F) is the differential operator defined in Definition 13.61 and {C; F) is the Berline- 
Vergne operator defined in Definition II. 1[ 

Proof. For any rational space (X, A), any rational subspace L in X, any unimodular cone C 
in X/L and any non-negative integer n satisfying n > dim(C), define the operator T>^{C) on X 
by the lift of Vn^^ (C ; 0) to X under the identification X/ L = L-^Q . We need to check that these 
operators satisfy the conditions in Theorem 11.21 The condition (1) in Theorem 11.21 follows from 
this definition. The condition (4) in Theorem 11.21 follows from Theorem 13.71 and Lemma 13.81 The 
condition (3) follows from Example above. The condition (2) follows from Definition 13.61 Therefore, 
the assertion follows from Theorem 11.21 □ 
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4 Asymptotic Euler-Maclaurin formula over rational cones 

In this section, we derive an asymptotic Euler-Maclaurin formula of Rn{C](p) for general rational 
cone C. To discuss asymptotic expansion of R^^C^cp) for pointed rational cones C, we define, for 
such a cone C and non-negative integer n, the distribution An{C; •) G <S'{X) by 

An{C;ip):= Yl / D^(C;F)ip, (4.1) 

FGJF'(C),dim(F)>dim(C)-n ^ 

where D^{C;F) is the Berline-Vergne operator defined in Definition ll.il 

Lemma 4.1 Let {Ci}f^i be a family of pointed rational cones in a rational space (X, A) satisfying 
'^iTixiCi) = 0, where, for each subset S C X, xi^) denotes the characteristic function of S. We 
set m = maxjdim(Cj). Suppse further that there exists a vector r] G X* such that {r},x) < for 
each / a; G UjCj. Then, for each (p G S{X), we have 

^ ?-i^n-m+dim(C,)(Ci;V5) = 0. (4.2) 

i,dim(Ci)>?Ti— n 

Proof. Since Ak{Ci; •) are distributions and Cq^{X) is dense in S{X), it is enough to prove ()4.2p 
for each ip G C^{X). Note that the function S{C) defined in (|l.ip have a valuation property (see 
[bp]). By this and the equation (11. 2p . we have 

E r,MvrG(Q))/(G)=0, 

where the subscript X/L{G) in iix/l{g) is dropped since these functions are lift to X* . Substituting 
(t G M, ^ G X*) in these functions and taking the Taylor expansion of each function, we have 

Y E t^^i^'^^''^'-H^G{Ci)mI{Gm = 0. 

k>-m i,Ge J"(C,),dim(G)+fc>0 

Thus, each coefficient of t^ in the above vanishes, and hence we have 

Y r,;u"-™+^^-(^)(^G(C.))(iC + il)I{Gm + r/) = (4.3) 

i,G€r{Ci),dira{G)>m~n 

for each n > 0, where, r] G X* is as in the statement of the lemma and ^ G X* is arbitrary. Let 
93 G C^{X). We have 

Dn-^+6iHc.){Cr;GMx) = ^-i^^^e,5+,(x)/."-'"+'i^-(«)(7rG(C0)(ie^ (4.4) 

where the Lebesgue measure d^ on X* is normalized as in Subsection 11.31 Taking the integral over 
G, we have 

/ I?n-,n+dim(C.)(Q; G)^ = (2^)— / I{G){ii + r?)^"-'"^'^''"^^) (^G (Q)) (i^ + " ^C, 

(4.5) 

where we have used the fact that ej^+^ is integrable on G for each i and G G T{Ci). Thus, 
multiplying (j4.5p by r^, taking the sum over all i and G G T{Ci) with dim(G) > m — n and using 
the equation (jl3D, we have (jOj) . □ 
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Theorem 4.2 Let C be a pointed rational cone in a rational space A) with a rational inner 
product Q. Then, for any ip € S{X), we have 

n>0 

where An{C;f) is defined in (j4.ip . Furthermore, the uniqueness statement of Theorem 11.21 still 
true if we replace the unimodular cones in the statement of Theorem 11.21 with the pointed rational 
cones. 

Proof. By replacing X with L{C), we may assmne that m := dim(C) = dim(X). It is weh-known 
that, for any pointed rational cone C, one can find a finite set of unimodular cones C = {ai}f^^ 
such that C is a subdivision of the pointed cone C, namely, the collection C satisfies the following. 

(1) C = y -F(cj) (2) (T, T G C ^ a n T G J-(cj) n -F(r) (3) C = \J a 

o-eC o-eC 

(For a proof of this fact, see [FuJ, Section 2.6.) By the inclusion-exclusion principle, there is a 
relation x{C) = J2aeC ^o-xi^^) with some G Z. Then, we have 

o-eC 7eA aec 

Each cone o" G C is a unimodular cone in X, and hence we can apply Theorem 13.71 to i?Ar(o"; (/?) for 
each a £ C. By a direct computation, we have 

n>0 (tGC ,dim(a)>m—n 

and hence Lemma 14.11 shows the first part of the assertion. The last assertion on the uniqueness 
follows from the same discussion as in Theorem 11.21 and hence we omit the proof. □ 

In the next section, we need the following lemma, which generalizes Lemma 13.21 

Lemma 4.3 Let C be a rational cone in a rational space (X, A). Let L = C D (— C). Then, for 
any ip G S{X), we have 

Rn{C; ^) = RNiTTLiC); {ttl^v) + 0(A^-~). 

Proof. If L = {0}, we have the conclusion without the term 0{N^°°). So, we assume that 
L ^ {0}. For simplicity, we write vr = vr^ : A" — > X/L, the natural projection. Take ip G Cq°{X). 
Since L is rational, one can take a complementary rational subspace W to L such that X = L(BW 
and A = (L n A) © {W D A). Set G = C nW , which is a pointed rational cone in W. We have 
C = L + G. Take a subdivision C of G into unimodular cone in W . The set {G^ = L + a ; cr G C} is 
a subdivision of G into rational cones. Then, there is a relation x(C) = So-eC '^o-x(Co-), and hence 

RN{G;ip) = J^iV— +'i^-(^'^)r.i?^(a;<^). 

a-ec 

Note that the cones G^ is of the form discussed in Lemma [3^ By Lemma [3^21 we have Rn{Go-] ^p) = 
RNiiT{cr);'iT^ip) + 0{N-°°), and hence 

RN{G;ip) = ^iV-™+^*"^('^)+'i''"(^V,i?^(^(a);^,(^) + 0(iV-°°). 
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The set {7r(cr) ; cr € C} is a subdivision of the pointed rational cone vr(C) in X/ L into unimodular 
cones. Furthermore, since xiP) = X^o-eC ^^"/^(^o") have x(^(C)) = X^o-eC ('^ defines 

a valuation. See [BP].) Thus, the sum in the right hand side of the last equation is i?jv(7r(C); tt^kC/?), 
which proves the assertion. □ 



5 Results and their proofs 

In this section, we restate Theorem [1] on the asymptotic Euler-Maclaurin formula of the Riemann 
sum 

7e{Arp)nA 

for a lattice polytope P in a rational space (X, A) and a smooth function on P in the abstract 
notation we used as before and give its proof. We also state and give proofs of its corollaries. 

5.1 Main theorems and their proofs 

Theorem 5.1 Let P he a lattice polytope in a rational space (X, A) with a rational inner product. 
For each f G J^{P) and n € Z+ satisfying dim(/) > dim(P) — n, let [P] f) be the differential 
operator defined in Definition 11.11 Then, for each ip € C°°(P), we have the following asymptotic 
expansion: 

P^(P;(^) ~ J]A„(P;^)7V-", A„(P;^)= /"z^^(P;/)(^. (5.1) 

n>0 /GJ="{P);dim{/)>dim(P)~n •'^ 

To prove Theorem 15.11 we need the following lemma. 

Lemma 5.2 Let f G -P(P) and let n G Z+ satisfy dim(/) > dim(P) — n. Then, for any g E J~{P) 
such that g C f , we have 

D^{P-f) = D^{^g{Cp{g))-7rg{Cf{g))). (5.2) 

Proof First of all, as in Subsection O note that we have D^{P;f) = {TTf{Cp{f));Q). We 
set C = 'Kg{Cp{g)) and G = iTg{Cf{g)). Then, C is a pointed rational cone in X/L{g) and G € 
F{G). Furthermore, we have {'Kg{Gp{g));'Kg{G f{g))) = {t^g{G);Q), where tig ■ X/L{g) -> 
{X / L{g)) / L{G) = X/L{f) is the natural projection. Since 'KG°'^g = '^f '■ ^ ^ ^/L{f) and 
Cp{f) = L{f) + Gp{g), we have 7rG(C) = ■Kf{Gp{g)) = tt f{Gp{f)), and hence the equation (j5.2p 
follows. □ 

Proof of Theorem 15.11 For any g G -P(P) and v G g, we set Gp{g) = Gp{g) + v which does 
not depend on the choice of u G 5. Then, we use the following version of Euler's formula ( |BrVj . 
Proposition 3.2, (1)): 

5{{NP)r^K)= J2 {-lf'^^''^S{C^p{Ng)nA), (5.3) 
seJ-{P) 

where, is a positive integer and, for any subset S of A, 5{S) is a distribution defined by 
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For each N £ Z>o and G C°°(X), we set {D*yj^ip){x) = ^{x/N). For each g G J^{P), we fix 
Vg €z g n A. Clearly we have 

( 6{NP n A),Dl/^^ ) = iV^i-(^)i?;v(P; 
{6{C+p{Nf)nA), Dl/^^) = N'''^^P^RN{Cp{g)-T,^^), 

where, for f G X, we set Tyip{x) = ip{v + x). Take if G C°^{P) and extend 93 as a compactly 
supported smooth function on X. Then, by (j5.3p and Lemma we have 

^'^^"^ . (5.4) 

^ ^ (-l)'i-(^)i?^(vr,(Cp(5));(vr,).r,,V^). 

Since ■Kg{Cp{g)) is a pointed rational cone in X/L{g) with respect to the lattice TTg{A), we can use 
Theorem 14.21 for i?jv(7rg (Cp(g)); ^-^^d hence 



dim ( G ) > dim ( P ) — 71 — dim (3) 



n>0 

(5.5) 

G 



MP;^)= E i-lf^'-^^ [ D^{7:g{Cp{g));G){7rg),T, 



Each faces G G F{-Kg{Cp[g))) with dim(G) > dim(P) — n — dim((jr) can be written as G = 7rg{Cf{g)) 
with a face / G J~{P) such that g C f and dim(/) > dim(P) — n. Furthermore, the correspondence 

{/ G J-(P) ; 5 C /} 9 / ^ %(C/(<7)) G J-(%(Cp(5))) 

defines a bijective correspondence between the above two sets. Thus, by Lemma 15.21 and the 
definition of the function (TTg)^T^j,ip, we can write 



A(P;(^)= Yl E i-ir^^'-'^ [ D^{P;f){ng).n^^ 

dim(/)>dim(P)— n 

= E E i-ir^'-^^ D^{P;f)^ 

dim(/)>dim(P)-ra 

tt;-. . ~;:7,s J(f) 



f€HP), dim(/)>dim{P)-n geT{f) 

where ( / ) is the affine hull of /, and for each S C ( / ), we denote X'^i'S) the characteristic function 
of S on ( / ). In the first line above, we used an obvious identity D^{P; /)(%)* = {'^g)*Dn f )"^ 
for ip G C^{X). To simplify the above, we use the formula (Proposition 3.1, (1) in [BrV] ) 

E i-^f''''^'^x{C+{g))=x{P). (5.6) 
seJ-{P) 

Note that in [BrVj . the above formula is proved for P with non-empty interior. Replacing P by 
/ G J^{P), which is regarded as a polytope in the affince subspace (/) with non-empty relative 
interior, we have 

E (-i)''"^'^x^(c^/(5)) = x^f). 

g&Hf) 
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Therefore, we obtain the formula (|5.1|) for An{P; which complete the proof of Theorem 15. II □ 

Remark: In the proof above, we used Theorem 14.21 However, if the lattice polytope P is Delzant, 
then the cone TTf{Cp{f)) for each / € J'{P) is a unimodular cone in X/L{f). Therefore, we only 
need to use Theorem 13. 7[ Hence, for Delzant lattice polytopes, it turns out that our proof of 
Theorem 15.11 is independent of |BeVj . However, for general lattice polytopes, it does not seem to 
be easy to construct the operator D^{P; /) in such a way given in Definition l3.6l Indeed, Definition 
13.61 is based on Proposition 13.11 This means that if we could obtain a result like Proposition 13.11 
for general rational cones, then one might be able to find such an expression as in Definition 13.61 
Hence, it might be better to prove a result like Proposition 13.11 for rational cones without using a 
subdivision of rational cones into unimodular cones. However, to do this, it seems that one need 
to find a different method. 

Next, we show that, under some assumptions, the asymptotic expansion of Rn{P',<p) of the 
form (jS.ip is unique. 

Theorem 5.3 Suppose that, for any rational space {X, A) with a rational inner product, rational 
subspace L of X, pointed rational cone C in X/L and non-negative integer n such that n > dim(C), 
there exists a homogeneous differential operator 'D^{C) of order n — dim(C) with symbol u^{C) 
such that they satisfy the conditions (1), (2) and (3) in Theorem 11.21 Furthermore, suppose that, 
for any lattice polytope P in X and ip € C°°{P), the following holds: 

Rn{P; ^) ~ iV-" Yl I T^nMCpif)))^. (5.7) 

«>0 /eJP'(P);dim(/)>dim(P)-n •'^ 

Then, we have D^iC) = D^{C;0) for any pointed rational cone C in X and non-negative integer 
n with n > dim(C), where the operator D^{C;0) is defined in Definition ll.il 

Proof Let us prove the assertion by the induction on dim(X). For dim(X) = 0, 1, the assertion 
is true by the condition (3) in Theorem 11.21 Suppose that for each {X,A) with dim(X) < m — 1, 
the assertion holds. Let dim(X) = m. Take a pointed rational cone C in X. We may assume that 
dim(C) = m. Take a vector ^ G A* such that {£,,x) > for any x € C. Set Pi = C R {x ; {£,,x) < 
1}, which is a rational polytope in X. Hence, each vertex of Pi is a rational point in X. We 
take a positive integer q such that P = qPi is a lattice polytope. Let U he a small open ball 
around the origin such that U n V(P) = {0} and U C {x ; {S.,x) < q}. Then, by the assumption, 
for each (p € C^{U), the Riemann sum R]\f{P;(p) admits the asymptotic expansion (15. 7p . In 
()5.7p . if dim(/) > 0, then since dim(7rj(Cp(/))) = m — dim(/) < m, the differential operators 
V^{'Kf{Cp{f))) coincide with D^{PJ) = D^{-Kf{Cp{f)),0) by the induction hypothesis. Take a 
vertex v of P. Suppose v ^ 0. Since is zero near v, the contribution from the vertex v to the 
expansion (15. 7p vanishes. Thus, by Theorem l5.lt we have 

[V^{MCpm)vm = [D^{P;0)vm in > m). (5.8) 

Take p G C^{U) such that p = I near 0. For any ip G C~(X), we have I^M) for ip = pip. But 
since p = I near 0, the equation (|5.8p holds for any ip G C°°{X). Take ip G C°°{X) and x G X. 
Applying (15. Sp for the function T^ip, we have 

[D^{no{Cp{0)Mix) = [D^{P;0)ip]{x) {n > m) 
for any ip G C°°{X). Since 7ro(Cp(0)) = C, we conclude the assertion. □ 
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5.2 Computation in one and two dimensions 

A polytope P in a rational space {X, A) is said to be Delzant if for each vertex v of P, the number 
of edges incident to v is dim(X) and there exists an integral basis ii^ of A such that each edge 
incident to v is of the form {v + te; t > 0} with an e € i?. In this and the next subsection, we give 
explicit computations for Delzant lattice polytopes. To compute each coefficient A„(P;(/9) in the 
asymptotic expansion of the Riemann sum Ri\f{P; ip), it is important to compute in low dimensions. 
In this subsection, we perform these computation. In this and the next subsections, we drop the 
superscript X in D^{P,g). 



5.2.1 In one dimension 

Let X he a. 1-dimensional vector space with the lattice A. Let n € A be a generator and set 
C = M+ii. We have computed the differential operator Dn{C; 0) in Example at the end of Subsection 
13.31 Let P be an interval given by P = {tu G X ; a < t < b} with a,b a < b. Since 

Dn{P; P) = for n > 1, we have 

An{P; = DniP; {a}Ma) + Z)„(P; {b}Mb) = (-l)"-ip(n)[Vr V(a) + V!^;V(6)] 
Identifying X = M and n = 1 so that A = Z, we have 

n! 

Substituting b2m+i = (m > 1) and 62m = (— l)™'~^Pm with the Bernoulli number Bm, we have 

[2m)\ 

which shows the classical asymptotic Euler-Maclaurin formula. 



5.2.2 In two dimension 

Next, we compute in two dimension. Let (X, A) be a two dimensional rational vector space with 
a rational inner product Q. Let E = {61,62} be an integral basis of the lattice A, and set C = 
M+ei +M+e2. Set 

61 = til + 6162, 62 = ^2 + 6261, 

where the non-zero vectors ui,U2 € X satisfy Q(ni,62) = Q{u2,ei) = 0, and 61,62 E Q are given 
by 

Q(ei,e2) Q(ei,62) , . 

ci = r, C2 = r. (5.9) 

Q{e2,e2) Q{ei,ei) 

Define Ai, A2 € by Xi{ej) = 6ij. A straightforward computation shows 

L{C;E,E;kXi+lX2) = Vty,^, 
L{C- {ei}, {ei}; fcAi) = V^;, , L{C- {62}, {62}; /A2) = VL,, 

k-l l-l 

L{C- {ei}, E- kX,) = c^Y, V^, V,Y'-^ L{C; {62}, E; IX2) = 62 ^ V^^^-^ 
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Set Fi = M+e2, F2 = M+ei. Then, we have 

Dn{C; Fi) ={-ir-^p{n)Vl;\ Z)„(C; F^) = {-ir-Mn)K;' (n > 1), 

n-l 

Dn{C; 0) =(-1)" 5]p(A;)p(n - fe)V,Y'V^-i-^ 

k=l 

(n-2 n-2 \ 

ci E +C2Y ^lyV'' > 2). 
s=0 s=0 / 

Let P be a Delzant lattice polytope in (X, A). For each facet / of P, Dn(P; f) is the hft of 
Dn{'Kf{Cp{f))]Q). Let a/ G A be the inward primitive normal of /. (Such a vector aj exists 
because the dual basis of an integral basis of A with respect to Q is rational.) We identify TTf{Cp{f)) 
with R+aj by the map 

Let ei G A be a generator of L{f) fi A. Since P is Delzant, we can find 62 S Cp{f) D A such that 
{61,62} forms an integral basis of A. Then, the vector 

is a generator of (pf{Trf{A)) such that ipf{7Tf{Cp{f))) = M+uj. Note that the definition of ttj does 
not depend on the choice of 62 G Cp{f) fl A whenever 61,62 forms an integral basis of A. Hence, 
by (j3.26p . the differential operator Dn{P; f) is given by 

D^{P;f) = {-ir-Mn)V:;' = -^Kj' {n > 1). 
Therefore, we have the following. 

Corollary 5.4 Let {X,A) be a two dimensional rational vector space with a rational inner product 
Q. Let P he a Delzant lattice polytope in (X, A). Then, the coefficients An(P;{p) (n > 2) in the 
asymptotic expansion (j5.ip of the Riemann sum Rn{P; ip) is given by 

AniP;ip)= Y [Dn{P;f)^+ Y DniP;vMv). 

In the above, the differential operators Dn{P\f) and Dn{P;v) are given by 
r, (p. f \ — _^v""^ 



n-l 



bkK-k V7fc-1 wn-l-fc 



r) (p..,\-Y^ "k^n-k ^k-l y 



k= 

n-2 n-2 



\ s=0 s=0 / 



where, for a face f G F{P)i, G denotes the inward normal defined in ()5.10p . and for a 
vertex v G V{P), the vectors ei{v),e2{v) G A denote the integral basis of A such that two facets 
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meeting at v lie on the half lines v + tei{v), t > 0, i = 1,2, and ui{v),U2{v) € X satisfy 

Q{ei{v),e2{v)) 



ei{v) = ui{v) + ci{v)e2, Q{ui{v),e2{v)) = 0, ci{v) 
e2iv) = U2{v) + C2{v)ei, Q{u2{v),ei{v)) = 0, C2{v) 



Q{e2{v),e2{v)) 
Q{ei{v),e2{v)) 
Q{ei{v),ei{v))' 



Note that, in the following, we use D2{C;0) for two dimensional unimodular cone C. The explicit 
formula for D2{C; 0) is given by 

D2{C; 0) = p{lf + (ci + C2)p{2) = i + (ci + ca)^, (5.11) 



where ci , C2 are given in (|5.9p . 

5.3 Computation of the coefficient in the third term 

Our main Theorem 15. H or rather the construction of the operators Z?„(P; /), allows us to compute 
the coefficient A2{P; f) in the third term of the asymptotic expansion (j5.ip . Before computing the 
third term, let us compute the first and second terms. 

Corollary 5.5 For any Delzant lattice polytope P in a rational space A) with a rational inner 
product Q, we have 



where the integration on facets g € F{P)m-i is performed with respect to the measure on g induced 
by the lattice A. 

Proof. The first term is obvious. For the second term Ai(P; (p), note that the dimension of faces 
which contribute to Ai{P;Lp) is m — 1 and m. But the operator D\[P;P) is the lift of L'i(0;0) 
(see Definition 13. 6p which is zero. Thus, the contribution to Ai (P; ip) comes from facets. Let 
g G T{P)m-i- Then the operator Di{P;g) is the lift of Di{TTg{Cp{g));0), which is a rational 
constant. Let G A be inward primitive normal of the facet g. As in the computation in two 
dimension, let tpg : X/L{g) — )■ L(g)-^Q be the isomorphism defined by 



We take an integral basis ei, . . . ,em-i of L{g) n A. Since P is Delzant, one can take G Cp{g) 
such that ei, . . . , em form an integral basis of A. We set 



have Di{TTg{Cp{g));0) = — ^ = |. Hence, we have 



As before, the definition of Ug above does not depend on the choice of Cm above. By (j3.26p . we 

lence, we have 



m — 1 
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which completes the proof. □ 

Note that the above formula for the second term Ai{P; ip) coincides with that in (j0.7p . Indeed, 
if X = M'", A = Z"* and Q is the standard Euclidean inner product, the primitive inward primitive 
normal Ug for each facet g oi a Delzant polytope P is a part of an integral basis of Z™". 

Next, we compute the third term, which does not seem to have been obtained before. For 
simplicity, we work in the Euclidean space X = R*" with the standard lattice and the standard 
inner product. 

Corollary 5.6 Let P he a Delzant lattice polytope in the Euclidean space (]R'",Z"') with the stan- 
dard inner product Q. Then, we have the following: 



1 1 f Q{ai{g),a2{g)) ^Q{ai{g),a2{g)) 



4 12 \Q{ai{g),ai{g)) Q(a2(ff), 02(5)) 



where, for g E F{P)m-i, the vector Og is the inward primitive normal to g, and for g E T{P)m-2, 
the vectors ai{g), a2ig) are the inward primitive normal to the facets (71,52 £ ^{P)m-i such that 
9 = 91 n52- 

Proof. By (|5.ip . the faces which contribute to A2{P; 99) is of m — 1 or m — 2 dimension. Let g 
be a facet of P. Then, Dn{P;g) is the lift of Dn{TTg{Cp{g));0). Hence, as before, we have 



D4P;g) = {-ir~Mn)K;' = '^-^"^V ^ 1) 



where the rational vector Ug G L{g)^Q is given in (j5.12p . But, we are working in the standard 
Euclidean space with the integral lattice Z™" and the standard inner product. Since P is Delzant, we 
can take an integral basis ei, . . . , of Z™ such that ei, . . . , Cm-i is an integral basis of L{g) n Z™ 
and if we denote the dual basis of ei,...,em by qi,...,q^, then = cxg. Thus, we have 
Ug = ag/Q{ag,ag) and hence 

Next, suppose that g is a face of dimension m — 2. Take two facets gi, g2 such that g = gi H g2- 
Denote ai{g) G A the primitive inward normal to gi {i = 1,2). Let f be a vertex in and take 
gs, ■ ■ . ,gm £ J^{P)m-i such that {v} = 51 n • • • n gm- Let E = {ei, . . . , Cm} be an integral basis of 
Z™" such that each vector v + ej defines an edge incident to v and v + ej ^ gj . We have 

Cp{g)=R+ei+R+e2 + L{g), 

and 63, . . . , Cm is an integral basis of L{g) n Z™'. Let ai, . . . , Um be the dual basis of ei, . . . , em- 
Then Oi = Oii{g) for i = 1, 2, and ai,a2 form a basis of L{g)^ . We write 

ei=tii+fi, e2=U2 + V2, ui,U2&L{g)^, vi,V2 £ L{g). 

Under the identification 

X/L{g) 3 x + L{g) ^ Q{x, ai)ui + Q{x, a2)u2 G L{g)^ , 
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the cone Trg{Cp{g)) is identified with M+ui +M_|_U2 and the generator of 7rg(Z"^) is identified with 
ui,U2- Thus, by (|5.1ip . we have 

D2{P;g) = D2{TTg{Cp{g));0) = - + — — + — . 5.13 

4 12 \Q{ui,ui) Q[U2,U2)J 

But then it is straight forward to show that 

Qia2,a2) . Q{ai,ai) Q{ai,a2) 
Q{ui,Ui) = — , Q[U2,U2) = — , Q{ui,U2) = — , 

D = Q{ui,Ui)Q{u2,U2) - Q{ui,U2f- 

From this and (j5.13p . we conclude the assertion. □ 
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